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In [l] we suggested the factorization formula for the nucleon form factors which consist of the 
sum of two contributions describing the hard and soft soft spectator scattering, and we provided a 
description of the soft rescattering contribution for the FF _Fi in terms of convolution integrals of 
the hard and hard-coUinear coefficient functions with the appropriate soft matrix elements. 

In present paper we investigate the soft spectator scattering contribution for the FF Fi. We 
focus our attention on factorization of the hard-coUinear scale ~ QA corresponding to transition 
from SCET-I to SCET-II. We compute the leading order jet functions and find that the convolution 
integrals over the soft fractions are logarithmically divergent. This divergency is the consequence of 
the boost invariance and does not depend on the model of the soft correlation function describing 
the soft spectator quarks. Using as example a two-loop diagram we demonstrated that such a 
divergency corresponds to the overlap of the soft and coUinear regions. As a result one obtains large 
logarithm InQ/A which must be included in the correct factorization formalism. 

We conclude that a consistent description of the factorization for _Fi implies the end-point collinear 
divergencies in the hard and soft spectator contributions, i.e. convolution integrals with respect to 
collinear fractions are not well-defined. Such scenario can only be realized when the twist-3 nucleon 
distribution amplitude has non- vanishing end-point behavior at a certain low normalization point. 
Such behavior leads to the violation of the collinear factorization for the hard spectator scattering 
contribution. 

In order to perform factorization of the hard-coUinear modes one has to formulate an unambiguous 
method for the separation of the soft and collinear modes in SCET. We expect that the soft spectator 
scattering plays important role in such derivation. 

We also discuss the physical subtraction scheme for SCET-I factorization which can be used for 
systematical analysis of the hadronic processes in the range of moderate values of ~ 5 — 20 GeV'^ 
where the hard collinear scale ~ QA is still not large. 



I. INTRODUCTION 

A substantial progress of the experimental studies of the nucleon form factors (FFs) has been achieved 
during the last decade. The polarization transfer method allowed one to measure accurately the proton 
FFs up to momentum transfer ~ 8.5 GeV^ [2HS] ) for recent reviews see, e.g., Refs. [SHE]. It also 
opened a possibility for systematic studies of the FFs at large space-like region in the near future at 
the Jlab 12GeV upgraded facility see, e.g., [5]. At the same time the PANDA Collaboration at GSI is 
planning to carry out precise measurements of the proton FFs at large time-like momentum transfers, 
up to around 20 GeV^, using the annihilation process p + p e+ + e~ [TOl E]. These experiments will 
provide us with new information on the FF behavior in the region of large momentum transfers that also 
provides a strong motivation for the theoretical studies of the large-Q behavior of the nucleon FFs. 

It is known for a long time that in the case of nucleon FFs the soft spectator scattering^ is not 
suppressed by inverse powers of large Q [T^HTS] . 



* On leave of absence from Petersburg Nuclear Physics Institute, Gatchina, 188350, Russia 

^ Let us mention that in the literature there are many equivalent names associated with the soft spectator scattering: 
double flow contribution (or regime), soft-overlap contribution, Feynman mechanism. 
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Moreover, results of different phenomenological approaches P^H22] allows one to conclude that in the 
region of moderate ~ 5 — 15GeV^ such mechanism plays the crucial role in order to obtain correct 
description of the FFs data and also other hard processes with baryons see, e.g., P5ti25j . A systematic 
QCD approach for description of such soft mechanism is needed in order to develop a model-independent 
formalism for description of existing and future experiments. 

In [T] the contribution of the soft spectator scattering has been included into factorization scheme 
using the formalism of the soft collinear effective theory fSCET) [2Sl - l5T| . The important feature in the 
description of the soft-overlap contribution is the presence of the hard-collinear scale ^ QA. In order to 
take into account the soft rescattering mechanism one uses the two-step matching technique developed 
in SCET. Following from the leading-logarithmic analysis we suggested that the tentative factorization 
formula for the Dirac FF Fi can be written as a sum of the hard and soft rescattering contributions: 

F, F^''^ + F^'\ (1) 

where the hard rescattering part F^'''^ is described by convolution of the hard coefficient function H with 
the nucleon distribution amplitudes (DAs) 

F^"'^ = J Dx, J Dy, ^{y,) H(x„y,;|g) = * * H * (2) 

while the soft contribution F^*-*-* has the same scaling behavior l/Q^ as the hard spectator term F^''-' 
and can be presented in the following form : 

Dy,^{yi) J dwidwa J'(y^,w^Q) J Dx,^{x,) J dvxdv^ 3{x,,v,Q)S{uj„v,). (3) 

This formula can be interpreted graphically as reduced diagram in Figjl] The hard subprocesses in 
the soft spectator contribution are described by the hard coefficient function Ca and two hard-collinear 
jet functions J, J'. They describe the parton scattering with the hard and hard-collinear momenta, 
respectively. The non-perturbative DA and soft correlation function (SCF) S describe the long distance 
scattering of collinear and soft modes. The convolution integrals in Eq. ([s]) are performed with respect 
to the collinear fractions Xi and y^, and with respect to the soft spectator fractions w^, Ui ~ A. 




FIG. 1: Interpretation of the soft rescattering as a reduced diagram 

In [1] we discussed the matching of electro-magnetic current onto SCET-I operator, calculation hard co- 
efficient function Ca and resummation of Sudakov logarithms. The structure of the factorization formula 
([3]) was considered only qualitatively. In the current publication we present more detailed consideration of 
this contribution. We shall carefully consider the matching of SCET-I to SCET-II and check the validity 
of the Eq.([3]). 

In ^1, we assumed the existence of the all convolution integrals with respect to the collinear and soft 
fractions in ([3|. Such assumption is motivated by the following observation. First, one can easily observe 
that collinear convolution integrals in both contributions ([2]) and (|3| have a similar structure. Second, if 
one uses phenomenological models of the nucleon DAs existing in literature then the collinear integrals 
are well defined. After that we expected that the integral over the soft fractions must be also well defined. 
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However in this work we show using some general model-independent arguments that the convolution 
with respect to the soft fractions in ^ is divergent. This divergence can be represented as scaleless 
logarithmic integral over one of the fractions; ^ duji/uji. Obviously, such integral has problems in 
both ultraviolet (UV) and infrared (IR) regions. The divergence in UV-limit is a signal that here one 
faces with the well known problem: overlapping of the soft and collinear modes. 

In order to clarify the situation we perform more careful analysis of the factorization ([T]) in the per- 
turbation theory using perturbative expressions for the DAs and soft correlation function. It turns out 
that in the theory with massive quarks the collinear convolution integrals in ^ and ^ are also diver- 
gent. This divergency arises due to the overlap of the soft and collinear sectors. We demonstrate that 
there is specific large logarithm ^ InQ/A which can be computed due to the overlap of collinear and 
soft regions associated with ^ and (|3]).This is exactly the logarithm which was computed in [T^. In 
order to compute ([2| and ([3| unambiguously one has to define a certain prescription which allows one 
to separate the collinear and soft sectors and to avoid a double counting. In the case of perturbation 
theory such separation can be carried out using the the idea of subtractions discussed in QCD for the 
Sudakov FF [32] or, similar technique in SCET known as zero-bin subtraction method [33]. We obtain 
that the perturbative calculations of the studied diagram is in agreement with the factorization formula 
([2]) . However generalization of the perturbative results to realistic case faces with certain difhculties and 
the description of the hard-coUinear factorization even for Fi remains challenging. 

The overlap of the collinear and soft regions imposes also a qualitative restriction on the end-point 
behavior of the twist-3 nucleon DA. We expect that the end-point behavior of nucleon DA at low normal- 
ization differs from the behavior of the asymptotic DA "^asi^i) ~ X1X2X3. The asymptotic-like models of 
DA always vanish when one of the fractions Xi tends to zero. But DA at a low normalization point should 
not vanish in order to produce the end-point singularity in the convolution integrals. We show that such 
behavior can be obtained even in perturbation theory with massive quarks: for instance, if 2:2 — )■ and 
xi is fixed, then from the calculation of the Feynman diagrams ^{xi, 0, 1 — xi) 7^ 0. We expect that the 
nonperturbative DA could posses the similar behavior. 

Our presentation is organized as follows. In Sec. II we describe in details the soft rescattering contribu- 
tion F^^^ defined by Eq.(b|, present the analytic results for the leading-order jet functions and discuss the 
divergency of the soft convolution integrals. Sec. Ill is dedicated to the analysis of the overlap between 
the soft and collinear modes. Using two-loop QCD diagram we calculate the "soft-coUinear" logarithm 
and perform the interpretation of the obtained results in terms of hard F^'^^ and soft F^*-*-* contributions. 
Generalizing these observations we suggest that collinear integrals in the factorization formulas for F^^'^'^ 
and F^^"^ are also divergent and therefore the nucleon DA must have nonvanishing end-point behavior. 
In Sec. IV we discuss the application of the SCET factorization for the phenomenological analysis at 
moderate values of Q^. In Sec.V, we summarize our results. 



II. SOFT RESCATTERING CONTRIBUTION 



We begin our discussion from the description of the soft rescattering contribution. All notations are 
the same as in [T] and are briefly described in Appendix A. For simplicity we consider proton as the 
target nucleon. In order to describe F^^"^ we need to define two nonperturbative functions: proton DA 
^'(xi) and soft correlation function (SCF) S{uji,i^i). 

The proton DA is a well-known object defined as: 



(4) 



where the measure reads Dxi = dxidx2dx36{l - 
index "c" denotes in SCET-II the collinear fields 



Xi - X2 



x^), the multiindex a = aia2'^3: £^nd the 



Wc[x\£,l = P exp <^ ig / dt {n ■ A,){x + tn) 



-Uc{x) 



(5) 
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The function "^aixi) can be further represented as 



75 



N 



-N 



7 75^ 



N 



(6) 



where C is the charge conjugate matrix C : C^^C = 7J. The scalar functions V, A, T depend on the 
colhnear fractions Xi and from the factorization scale which is not shown for simplicity. Alternatively, 
these functions can be represented through the one twist-3 DA (p]\f{xi) [33]: 



Vix,,X2. X3) = Vil, 2, 3) = -{^Nih 2, 3) + ^jv(2, 1, 3) ), 



A(l,2,3) 



--(^jv(l,2,3)-v^Ar(2,l,3)), 



T(l,2,3) = -(^^(1, 3, 2) + ^^(2, 3,1)). 



(7) 
(8) 
(9) 



In many applications ip^^Xi) usually is approximated by a few polynomials in Xi with unknown coeffi- 
cients. The polynomials represent the eigcnfunctions of the evolution kernels and naturally arise from 
the solution of the evolution equation ( see, e. 51., |351 136] and references there in). For example, the one 
of most popular parametrizations for ipN{xi) can be represented as follows (34j : 



ifiNixi) ~ 120xia;2a;3/jv(l + cio(a;i - 2x2 + X3) + cii{xi - x^)), 



(10) 



The values of the nonpcrturbative constants /n , c,; are not important for our further discussion and we 
will not provide their numerical estimates. 

A parametrization such as in Eq.( 10 ) always fulfills one important property: it vanishes at the boundary 
where one of the fractions is zero: 



ipNixi,X2 = 0, 1 - xi) = ^pn{xi = 0,X2, 1- X2) = (Pn{xi, 1 - a;i,X3 = 0) = 0. 



(11) 



Such behavior is an important requirement which ensures the existence of the convolution integrals in the 
factorization formula ([2| because the hard coefficient functions, as a rule, have the end-point singularities 
^ 1/2^1,2- The end-point behavior ( |11[ ) can be associated with the corresponding behavior of the evolution 
kernel. Performing the expansion of the DA Lpi^ in terms of the eigenfunctions of the evolution kernel and 
neglecting by the higher order harmonics one always obtains the model which vanishes at the boundary 
( 11 ). However such approach is consistent only if one assumes an appropriate convergence of the conformal 
expansion at a given normalization point. 

The second non-perturbative input is the SCF S introduced in Eq.([3]). In SCET it is constructed from 
the soft quark and gluon fields q and Ag , respectively. The gluon fields enter only in the form of Wilson 
lines such as 



Sn{x) = Pexp I ig I 

L J — oc 



ds n ■ Aa{x + sn) 



(12) 



The definition of S given in ^ implies that it is a tensor with respect to Dirac indices. It is convenient 
to rewrite it in terms of scalar functions. This can be done with the help of Fierz transformation so that 
the result reads: 



S{u!i, Vi) = ^Ct/i (g) ^fiC Sv{uJi, Vi) + 

00 o 

The scalar functions Sv,a,t defined as 

dri2 



o 4 4 



with the operators 



drji 
2ti 



2tt 



dXi 
'2^ 



dX 
Ytj: 



\Ox{^^,X^)\0), (14) 



1 k' 



5 



where 



X}X=V,A,T 



(16) 



The flavor structure of the SCFs Sx is defined by the flavor structure of the proton: it can be de- 
scribed as either uu- or Mc?-combinations. Below we will show that the tensor component St does not 
contribute to Fi and hence we can conclude that the matrix S{uji, Vi) is presented by four scalar functions: 



■'A.v 



(iOi, Vi). Indeed, SCFs also depend on the factorization scale, which is not shown for simplicity. 



One has to keep in mind that the arguments uji and vi (soft fractions) are of order A. Obviously, these 
variables can be associated with the light-cone projections of the soft spectators momenta which are 
defined to be positive. For instance, the leading order calculation of the SCF S'^{uji, Vi) in perturbation 
theory yields (see details in Appendix B) 



3m2 
167r6 



(17) 



where m denotes the mass of soft quarks. Notice that Sv,a is proportional to the square of mass rri^ 
which arises from the numerators of the soft quark propagators and this fact, as we will see later, has an 
important consequence. 

The last elements in Eq.([3]) which one has to introduce are jet functions J and J'. These are hard- 
collinear coefficient functions which can be computed in pQCD if the hard-collinear scale QA is quite 
large. They appear in the matching of hard-collincar modes onto collinear and soft fields in SCET-II. Jet 
functions can be computed from the T-products which schematically can be written as [1]: 



iC 



(fi) 

SCET-I 



(") 

SCET-I 



(18) 



where the asterisks denote the convolutions with respect to collinear and soft fractions. The matrix 
elements of the collinear operators in Eq.( 18 ) yield twist-3 DAs Q describing the initial and final protons, 
the soft fields q are combined into SCF S. Because the collinear operators in the both equations in ( 18 ) 



are the same, the functions J and J' are also the same, therefore one needs to compute only one of 
it. Corresponding leading order diagrams are shown in Fig [2] The last two diagrams with the three 





FIG. 2: Leading order SCET diagrams required for the calculation of jet functions. The inner dashed and curly 
lines denote hard-collinear quarks and gluons, external dashed lines correspond to collinear quarks, fermion lines 
with crosses denote soft quarks. Black square denotes the vertex of the SCET-I operator. 



gluon vertex have vanishing color factors and therefore do not contribute. This is in agreement with 
the observation made in Ref. [14] . The soft-coUinear vertices in the diagrams is obtained from the 
subleading SCET La grangian 'C^^^ see details in [T]. Computing these diagrams we projected the Dirac 
indices according to our definitions of the DA ^ and SCF (15). We find that the contribution with 
St corresponding to the tensor projection in eq. ( 16 1 has trivial coefficient function. This result is not 



an accident and can be explained by the angular momentum conservation (quite similar to discussion in 
[37]). Therefore we expect that corresponding hard kernel is trivial to all orders in as- 
In order to present the results for the jet functions let us rewrite Eq.([3]) in the following form: 



f['\Q) = Ca{Q) {e„ A"(Q) +ed /{'(Q)}, 



(19) 



where eu,d denote quark charges. The hard coefficient function Ca{Q) in Eq.(19) includes all large 



logarithms so that SCET-I form factors /"''^ depend only from the hard-collinear scale QA and defined 



as 



(P'ICW'7'^ W^Cb)sCET 



= N{p')'^l1_N{p) mQ,iiHc 0, QA) EE /f(Q), 



(20) 
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where the index q describe the flavor of the hard-collinear field (active quark). Computing the diagrams in 
Figj2]we obtained analytical expressions for the leading order jet functions. Our results can be presented 
in the following way 

f^{Q) = I'i{L0i,LJ2) * Sf{LJuUJ2,iyi,1^2) * 11(1^1,1^2) - * Sf{LJu^2,iyu^^^ (21) 

ftiQ) = 4K,^2) * 5r(wi,Wl,i^l,i^2) * 4(i^l,J^2) - /y(c^l,W2) * S^'\LJi,i02,iyul^2) * /y(;^l,i^2), (22) 

where asterisk denotes the convohition integral with respect to the soft fractions, for instance 

J^(wi,W2) * S'™(a;i,a;i,:/i,i^2) = / rfi:Jidw2 '^2)<S'^"(wi, wi, z^i, 1^2)- (23) 



The proton DAs and integrations over the coUinear fractions enter into the functions /y ^ in Eqs. ( 21 
These functions read: 



221 



Iy{uji,uj2) = (V - A-2T){yi,y2,y3)*Jaiyi,y2,y3,^i,i^2) 
+ [A-V - 2T){yi,y2,y3) * (2/1, y2, ys, ^i, ^2), 



(24) 



11(^1,(^2) = {V - A-2T){yi,y2,y3)*Jb{yi,y2,y3,uJi,uj2) 
+ (A-V - 2T){yi,y2,y3) * J^(2/i, y2, 2/3, ^i, ^2), 

Iy{ui,uj2) = V{yi,y2,y3) * ^a(yi, 2/2, 2/3,^1,^2) - A{yi,y2,y3) * (2/1, 2/2, 2/3, ^^i, ^2), 

/^(wi,W2) = -A{yi,y2,y3) * (2/1, 2/2, 2/3, ^i, ^2), 
where the asterisk again denotes the convolution integral over coUinear fractions: 

^(2/1,2/2,2/3) * ^'^(2/1, 2/2, 2/3, •••) = / -02/^ ^(2/1, 2/2, 2/3)^^^(2/1, 2/2, 2/3, •■•)• 



The leading order hard-collinear jet-functions J^f in Eqs.(24|27) read {yt — 1 — yi): 

87r2 



and 



•^"(2/1,2/2,2/3; Wi,W2) = al{^.hc) 

•^"(2/1,2/2,2/3; i^i,^^2) = alinhc) 

•^0(2/1,2/2,2/3; wi,'^2) = ali^ihc) 

•^^'(2/1,2/2,2/3; '^i,'^2) = alinhc) 



Stt^ 
27 

Stt^ 
27 



1 1 


[UJI - 


1- W2)a;iW2 2/22/3 




1 1 


(uji - 


\-u!2)ujiu!2 y2y3 




1 1 


[uji - 


\- U!2) UJ1IJJ2 2/12/2 




1 1 


{UJI - 


\-oj2)ujioj2 mm 



1 

2/1 
1 

2/2 



1 

^2 
1 

2/1 



4(a;i2/i +0^22/2! 
2/1 (wi -ha;2) 



(25) 
(26) 
(27) 

(28) 

(29) 
(30) 

(31) 
(32) 



The argument of the QCD running coupling is defined by the hard-collinear scale: ^uc ^ QA. From results 
(21 )-(22 1 one can deduce that the convolution integrals in (24)-(27l with the proton DAs vanishing at the 



boundary as in Eq.( 11 ) are well defined. Therefore one can assume that in the absence of other dominant 



regions, the convolution integrals with respect to soft fractions in (21|-(22| should also be finite. One 
can expect that the SCFs are concentrated in the region where the soft fractions is of order A and fall 
quickly in the region where the soft fractions are much larger than A. 



However the assumption of convergence is not correct. Using the perturbative expression (17) instead 
of nonperturbativc function one can easily obtain, that any convolution integral with respect to the 
soft fractions is logarithmically divergent and proportional to ^ ln[fii/v / I^ir] where scale /ijyv and /i/i? 
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represent the UV and IR cut-ofF respectively. This is exactly the logarithmic contribution which was found 
in [13] and later studied in [T3] and called as "nonrenormalization group type logarithmic contribution" . 
Using the leading order expressions given in Eqs.(p9|)-(32) one can show that the similar situation is 

ic soft fractions uJi and Vi can be associated 
Therefore the boost 



also takes place for the nonperturbative SCFs. Recall that t 

with plus and minus projections of the soft momenta describing soft spectators 
invariance implies that SCF depends on the products LOiVj: 



(33) 



Using this observation one can easily obtain that the convolution integrals of jet-functions and SCF in 
Eqs.(21 )-(22 1 are divergent. Consider, as example, the integral from the Eq.(21): 



J 



Using substitutions 



V2 



V2/tOl, i02^ Uj'2L0l 



(34) 



(35) 



one obtains 



Jq 



duj'o 



dv[ / dv'2ll{u:^,u'2UJl)Il{v[li^i,v'2/i^l)Sf{v[,v'2,L0W^,J2l^'2). (36) 



Using homogeneity of the leading order jet functions in Eqs.(29l-(32) yields 



I\{u;^,u:'2LOi)I\{v'Juji,v'2lu^) ^ I\{l,uj'2)ll{v[,v'2) 



and one obtains 



J = 



duji 



Jo 



duji, 



du[ / d^.^/;i(i,c.^)j:i(z.;,;.^)5^'i(:.;,^.^,c.^j.;,c.^z.^). 



(37) 



(38) 



The first integral in the last Eq.(38l is divergent Jg°° ~ \d\^uv / ^'Ir\- The similar arguments can be 
used for the all contributions in Eqs.(21 1-(22). Hence we can conclude that the soft convolution integrals 



are divergent and this singularity is independent on the particular properties of SCF. We performed the 
complete calculation of the jet functions in order to be convinced that this divergency does not cancel in 
the sum of all diagrams. 

It is natural to assume that dimensionless integral indicates that the soft region overlaps 

with the coUinear one as it happens, for instance, in the classical case of Sudakov form factor . 
Then we may also expect the appearance of the IR-divergencies in the collinear convolution integrals 
which describe the hard rescattering term and referred often as "end-point" divergencies. However, the 
corresponding collinear integrals computed with the existing models of DA, like one in Eq.(10l, are well 



defined. In order to clarify this situation we suggest to investigate the factorization properties in pQCD 
using as example appropriate 2-loop diagram which obtain contributions associated with the soft and 
collinear sectors. From analysis of the factorization of this diagram we may find the missing element in 
the discussed factorization picture. 



III. SOFT-COLLINEAR OVERLAP AND SEPARATION OF THE HARD AND SOFT 
SPECTATOR SCATTERING CONTRIBUTIONS 

A. Overlap of the soft and collinear regions 

In order to perform the required analysis we consider the diagram shown in Fig|3] This diagram was 
discussed in [T] in order to show nontrivial contribution arising from the soft region. But we did not 
investigate the possibility of the overlap of the soft region with other dominant regions. We consider this 
more accurately now because it will help us to solve the problem with the soft divergent integral and 
understand the factorization for Fi. The contribution of this diagram into the nucleon FF Fi can be 
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\j 



q=p--P 



FIG. 3: Two-loop diagram which obtain contributions from the hard, colhnear and soft regions. 



written as 



(39) 



where DA ^ is associated with the blobs in Fig|3] In order to simplify our consideration we use simple 
observation which follow from Eq.Q: 



'^{y,)i^Q, = 0. 



(40) 



This allows us in the intermediate calculations to substitute instead of full DA '9 the large components 
of colhnear spinors: 



*a(y«) 



4 



u{pi) 



_ 4 



4 

d{pz 



dips 



Kilo jella. Kale 



(41) 
(42) 



For simplicity we always assume that the bottom fermion line in the diagram corresponds to d-quark. 
We always assume that color indices are properly contracted and don't show them explicitly. Hence, 
contribution of the diagram in FigjS] can be associated with the perturbative diagram in Fig|4] and 
corresponding analytical expression reads: 



1 

o 
c 


2 o 


P 


2 p-p-k-k^% 
o 


o 
6 
c 
o 

rt> ° 


kj+k-p-p^ 



FIG. 4: Two-loop diagram with momenta flow. 



C / (rki(Pk2- 



k\ — ni?] [k. 



^i7* (fci + 'Tt) l^^i ^2l\p' -P3-ki+ m)7"(fc2 + ■m)^l^{p - ps - h + m)"f^ ^2 
[{P -P3-kif ~ m'^] Hp' -P's^ kif ~ m'^] {ki +k2-p' + p'^f (p - Ps - h ~ ^2)^ 

7" {p' - fci - fc2 + m) 7'' [p-ki- k + rnj -f^ ^3 
[(p' - fci - fc2)^ - m2][(p- ki - k2f - m?] {pi - kif [ki - p'lf 



(43) 



The quark mass m is used as a soft scale in order to regularize IR-divergencies and to describe the soft 
contribution, see Eq.(17|. The factor C accumulates color structures and others factors according to 
Feynman rules. 
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Our task is to clarify the relation between the collinear and soft regions performing the interpretation 
of I?^ in terms of different contributions according to factorization formulas ^ and If the overlap 
of these sectors is not possible then we may be find the missing elements in the factorization description. 
Notice that diagram in Fig|4] does not have UV-divergent subgraphs and therefore we do not need to 
consider renormalization of QCD parameters like quark mass and running coupling. Taking into account 



the logarithmic structure of tha hard spectator term 
the answer can be schematically written as 



{h) 



we may expect that after 2-loop integration 



1 



[A^^ In^ QVm^ + B'^-L In gVm^ + C^] + ©(l/Q^). 



(44) 



where the leading power is obtained from the dimension reasons. In the absence of the divergent 
soft contribution we could expect that coefficients A and B can be interpreted in terms of convolutions 
of the LO and NLO hard coefficient functions with the two- and one-loop evolution kernels, respectively. 
However, the divergency of the soft convolution integral in ^ allows us to suggest that there is a large 
logarithm of Q which has a different interpretation. 

The exact calculation of the coefficients A, B, C in Eq.([44| is a difficult task involving 2-loop massive 
integrals and we arc not going to do this calculation here. We will focus our attention on the possible 
overlap of the soft and collinear contributions which can most probably provide the solution of the 
problem. We will use the technique, known as strategy of regions [Ml EF! which allows to find and 
interpret the contributions originated from the different regions in ( 43 ) . This formalism also allows one 



to establish wether the collinear contributions overlap with the soft sectors. 
Let us begin our discussion from the soft region where 



kip. 



A, 



(45) 



Corresponding contribution can be represented as 



dk+ 



dk 



1,2 



xixl [-Q{k^ + k^)Y [-Qk^ 



dki2 



ki 



k2 + m 



02/32 



[kl 



[kl 



(46) 



The subscripts Ca,J,S are associated with the appropriate contributions in Eq.(|3| so that expressions 
in the brackets can be associated with these quantities. 

Computing the transverse integrals inside the brackets associated with SCF S one obtains the answer 



given in Eq.(17l which is proportional to m^. Redefining the light-cone variables fc^2 the convolution 



integrals one can eliminate the mass m from the consideration. But such redefinition performed in the 
divergent integrals may lead to a mistake. We will see that calculation of the regularized soft integrals 
leads to a logarithmic dependence on the mass m. 

We find that the contributions of the remaining regions in Z?^ can be associated only with the hard 
rescattering term F^'^'^ For brevity we skip the discussion of all the dominant regions and pass 

directly to those which are relevant to our consideration. These regions are collinear regions, where both 
spectator momenta fci^2 collinear either to initial p or to final p' momentum. Because of the symmetry 
between in and out states in Eq.(43) it is enough to consider only one of them. Let us choose 



kl 



p 



Q,k+^AyQ,kA 



A. 



(47) 



Then one obtains that corresponding contribution has the leading power suppression Q ^ as in Eq.(44) 
and can be written as a convolution of hard [...]h and collinear V parts (see details in Appendix C): 



c—p 



C [(7.^ 



dk^ dk^ V(fcj 



7-' 6 



Ixi [-Q{k^ + k^)]' [-Qfcr] 



(48) 
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where 



{^ItHp' - - fci + m)7"(fc2 + ™)} {el 7" {f 



ki 



{p' - ki - k2f' - m? {p' - p'^ - kif - m? {ki + k2 - p' + p'^f [ki - p'^)^ 



(49) 



The colhnear part V can be clearly associated with the two-loop contribution to the evolution kernel of 
the DA. Recall, that light-cone components k± scale according to (47). Computing the integrals over 
k'^2 ill Eq. ( 49 ) we obtain that the integration region for the minus components are restricted Q < k^ < 
Q. Therefore these variables can be rescaled to the dimensionless quantities which, as a rule, used as 



colhnear fractions. The integrals over the transverse momenta in (49) are UV-divergent and these are 
the logarithmical divergencies associated with the evolution of DA. 



In order to see the overlap of the soft and coUinear regions we consider the soft limit ( 45 ) in the colhnear 
contribution D^^ , . In this case one finds that the collinear-soft limit reproduces the soft contribution 
in Eq.(46): (i'c-pOs = D^^. This allows us to conclude that coUinear and soft regions overlap. As a 



consequence, this also allows us to expect that the convolution integrals with respect to fc_ in (48) are 



singular at the end-point region. We will see this explicitly later computing the colhnear integral (48) in 

secnmi 



Let us emphasize one important point. Computing the soft integrals in Eq.(46l we can skip the terms 
with momenta ki 2 in the numerator of the integrand 



dki I dky 



(l)fta, (1) 



[fc| — m^] [kf — m?] 



(50) 



This is possible because only m? contribution yields the final result, see Appendix B. From this ob- 
servation we can conclude that the part of the colhnear contribution D'^^^, in Eq.(48) relevant for the 
soft-collinear overlap is UV-finite. Thus we obtain that the soft-coUinear overlap in tiie given example is 
not connected with the UV-divergent colhnear sector which is associated with the evolution of nucleon 
DA. The importance of the mass squared term in the discussion of the soft spectator scattering was also 
noted in |13) . 

Taking into account this fact we can conclude that the convolution integrals of the hard coefficient 
function H in the factorization formula ([2]) with 1- and 2-loop evolution kernels Vi_2 or with the models 
of DA with the vanishing end-point behavior (11) (for instance, like (10)) are free from the end-point 
singularities. But the convolution of LO hard coefficient function H in Eq.(48l with the finite m^-term 



originating from the coUinear subdiagram V is singular. Below we will demonstrate that this singularity 
provides a logarithmic contribution which in the sum with the UV-logarithm from the soft term ( 46 ) 
yields a contribution proportional to InQ/m. 

Such situation turns out to be different from the many practical cases when the end-point singularities 
appear to be "visible" only after the inspection of the end-point behavior of the evolution kernel, or 
equivalently, asymptotic DA. Such situation takes place, for instance, in the case of Pauli FF F2. In 
this case the incompleteness of the colhnear contribution can be associated with the double logarithmic 
contributions. But in case of Fi the colhnear contribution has weaker singularity which is associated 
with subleading logarithm in the 2-loop diagram. As a result the hard rescattering mechanism can not 
describe correctly the coefficient B in Eq.(44). For its correct interpretation one has to involve the soft 
rescattering contribution. 



In order to clarify the interpretation of the UV-finite contribution arising in V in Eq. ( 48 ) consider 
the perturbative analog of the proton DA, i.e. as we would substitute three quark states instead of 
protons. In such case the 3-quark DA is defined as pQCD matrix element and its perturbative expansion 
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schematically read^: 



*pt(x,mj^) = *o(:c)+"^ *l(x,/iF)+ f"^^") *2(x,mf)+ ... . (51) 



where the higher order coefficients 5'i>o(a:;, (J-f) depend logarithmically on the renormalization (or factor- 
ization) scale fip. Schematically this can be represented as: 

*i(2:„/ij.) = ln/^|,/m2 Vi *«'o + 1'io(a;,), (52) 
'i'2{xt,HF) = InVfM^ ">^i*Vi**o + lnMFM^ "^2**o + ^'2o(a;i) + ..., (53) 

where asterisk as usually denotes the collinear convolution integrals, Vi 2 denote now 1- or 2-loop evolution 
kernels and dots represent the other contributions associated with the renormalization of QCD coupling. 
If we put /xf = m then all large logarithms vanish and we obtain: 

*PT(2:„m) = ^'o(x,) + ^^*io(a:.)+ M^ ^2o{x^) + ... ■ (54) 



Obviously, Eq.(54) is understood as a perturbative expansion of the DA at the low energy scale m ^ 
Q. Performing transition from perturbative consideration to physical FF one substitutes instead of 
perturbative DA the realistic one: ^pT(2^i,"T') — ^ ^{xi, Ho), where it is natural to assume that ^0 is 
a certain soft scale of order A. Following these arguments we can perform the interpretation of the 



UV- finite term ~ in Eq.(49). This term provides the contribution to ^'20 (2;;) and therefore can 
be associated with the DA at low normalization /ig. If '^2o{xi) has nonvanishing end-point behavior 
^2o(a^i,0, 1 — xi) 7^ then only the corresponding collinear convolution integrals are singular. Let us 
emphasize again that such a behavior is closely related with the presence of soft spectator contribution. 

Extrapolating this argument beyond perturbation theory we can not exclude that a similar situation 
takes place for the realistic DA. The singularity arising in the soft spectator contribution can be considered 
as strong argument in support of such scenario. In this case the nucleon DA (p]s[{xi, fj,o) at a low 
normalization point /io ~ A has nonvanishing end-point behavior. Such a behavior is very important for 
the consistent description of the large logarithm appearing due to overlap of the hard and soft spectator 
scattering mechanisms. With such DA one immediately sees that the naive collinear factorization ([2| for 
FF Fi does not hold anymore. 

Finally let us note that a rather similar conclusion about the end-point behavior of the nucleon DA 
was also obtained from the non-perturbative calculations. In pp] the nucleon DA ip^^Zi) has been 
computed in the chiral quark-soliton model [H] using large Nc expansion. The DA was estimated 
at ^0 ^ 600MeV and it was found that it does not vanish at the end-point limit. However in [40| it is 
also noted that obtained results are valid only in the region where ZiN^ ~ 1. Therefore for the realistic 
value Nc — i one can not expect an accurate description of the region with the relatively small collinear 
fractions. The nucleon DA computed at leading order in 1/Nc can potentially obtain large corrections at 
the end-point region. 

B. Calculation of the large logarithmic term arising from the overlap of collinear and soft 

regions 

Having established the soft-coUinear overlap one faces with the problem of unambiguous separation of 
these regions in the factorization formulas ([2]) and (|3|. This is quite a complicate problem, especially for 
the processes involving composite particles like hadrons. Therefore it is useful as a first step to consider 
more simple examples with similar features. We will continue the discussion of diagram in Fig|4] and 
discuss the separation of the soft and collinear contributions in this particular case. 



In order to compute the contributions of the soft ( 46 1 and collinear ( 48 1 regions explicitly we need a 
specific regularization for the convolution integrals over longitudinal momenta. Dimensional regulariza- 
tion (DR) can not be used in this situation because the soft contribution is given by scaleless integral and 



^ For brevity we skip the Dirac and color indices 
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therefore equals to zero in this case. Therefore one has to introduce a different regularization in order to 
work with the colhnear and soft integrals in ^ and ([3]). This is not only a technical problem. A careful 
prescription is required in order to avoid double counting computing the different contributions. These 
questions have been studied during last years in many publications in connection with the Sudakov form 
factor, see e.g. [321 |35] . In [32] a systematic subtraction procedure has been suggested in order to 
separate contributions of coUinear and soft modes in non-inclusive processes. The factorization in the 
context of SCET was discussed in [53] where the idea of the so-called zero-bin subtractions was invented. 
In our analysis we adopt this technique in order to formulate a prescription for the correct separation of 
the soft and colhnear modes and compute required integrals. 

Note that the soft contribution in Eq. ( 46 1 and the appropriate part of the colhnear term in Eq. ( 48 1 do 
not overlap with the hard region. As a result these integrals can be considered as UV-finite (in a sense that 
integrals over transverse momenta are finite). This allows us to carry out all calculations using a specific 
regularization in four dimensions which may be the simplest solution in this case. Such regularization 
prescription must be formulated uniformly for the all colhnear and soft divergent integrals. There is 
also a technical problem about applicability of the method of regions with this specific regularization in 
D ~ 4. In order to fix these details we suggest to investigate a simple one-loop integral which is close to 
our situation and can be easily computed. 



1. Colhnear and soft contributions in D = 4/ one-loop case 
As example consider following integral: 

/7T? 1 1 ?' 7r 

" Fr;;;3f iF^w^i = "'""^ + 

where we assume that expressions in the square brackets are defined with the +ie prescription: [X] = 
[X + ie], and the momenta p and p' are the same as used before. It is easily to see that this integral can 
be related to a well-known scalar vertex integral: 

" [fc2 _ m2] [fc2 _ 2{pk)] [fc2 - 2(p'fc)] ^^^^ 

The asymptote of the vertex integral is given by large Sudakov logarithm In^ /m^ but the mass differ- 
entiation reduces this structure to a simple logarithm. Analysis of the dominant regions yields: 



-'^coii ^ Js ^ g^, Jh ^ gj, (57) 



where the subscripts "coll", "s" and "h" denotes colhnear to p or p', soft and hard regions respectively. 
The hard region is suppressed and the large logarithm is generated only from the overlap of colhnear and 



soft regions. In order to obtain the leading order result (55) with the help of the method of regions we 
must introduce appropriate regularization. It is not difficult to see that DR can not help in this case. 

Consider the regularization by small off-shell external momenta introducing the small transverse com- 
ponents. This yields: 

[ m? I 1 

"^'■^^ ^ J'^^ [P-m2]2 [k^-2{pk)-pl] [k^-2{p'k)-p'l] ^^^^ 



As usually, one may expect that the exact answer ( [55[ ) can be reproduced by the sum of the contributions 
from the dominant regions: 

But in this case this not true and this can be easily seen from the explicit calculation. The contribution 
from the soft region reads: 

Js — I dk ^7 ; TT^-f ■ ^ — --r / dk 



[fc2 _ to2]2 [-p+k_ - pI] [-p'_k+ - pf] g2 J _ ^2]2 [_A:_ _ r_] - T+] 

(60) 
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where we introduced r_ = p\/p^ and = p''^/p'_. Notice that in the absence of the regulators this 
soft integral is scaleless and has both UV and IR divergencies as the integral in Eq. ( 38 1 . Performing 
integration over A;_ by residues and taking the transverse integral one obtains 



J. 



ITT 



dk: 



ITT m- 



In [tt'/t 



[k+ + t] [1 + k+T' /m? 



In [rr'/m^ 



(61) 



From Eq.(61 1 it is clearly seen that the regulators r± serves as IR- and UV-regulators {k+ — or fc_|_ — )■ cxd 
respectively). Passing to the last equation in (61 1 we neglected the regular in rr' contributions assuming 
tt' < vn?. 

In collinear region k ^ p' one obtains 



J c—p' 



= / dk 



1 



1 



dk^ 



m^f[-2{pk)-pl] [{p'-kY-p'l] 
dk I dfcj 



f-fc_ - • 



2] 2 [fc2-^2p'_fc+] 



(62) 
(63) 



We neglected the second regulator t+ in (63) because the corresponding integral is finite. Again, com- 
puting the integrals over k+ and fc_L wc obtain 



J c—p' 



it: 
P+ Jo 



[k- 



(64) 



Similarly one computes the second collinear integral k ^ p: 



J c—p 



(65) 



However the sum of the all terms (591 can not reproduce the exact answer (55 1. The reason is that 
collinear integrals Jc-p' and Jc~p obtain contribution also from the soft region. Computing the soft limit 
in expression in ( 63 ) 



1 



[/c2 



^^f{-2{pk)-p\\ {-2{p'k)-p'J\ 



(66) 



we reproduce expression for the soft integral in ( |61[ ). Hence in order to compute the contribution of 
the collinear regions correctly one has to subtract from the expressions for collinear integrals Jc-p' and 
Jc-p appropriate soft contributions [35] (or similarly to perform zero-bin subtractions [33] )■ With such 
subtractions one has 



\'J c—p' 'J s\ 



In 



'c—p 



ITT p+T^ 
- J., = in TT 



(67) 



Notice that logarithms in ( |67[ ) originates from the UV-region because IR singularities cancel in the 
difference exactly as it was discussed in [33] . Now the sum of the all contributions reproduces the correct 
answer: 



J — [<-^c— p' J s\ \Jc—p 'J s\ ^" Js ~Py2 2 



m J P+T- 



— In^ 

771,2 



(68) 



One can easily see that zero-bin subtractions ( 66 1 can be also taken into account by changing the sign of 
the soft contribution in Eq. ( 59 1 that was noted already in [33] • 

From considered example we can conclude, that evaluations of integrals using the method of regions 
in D = 4 with specific regularization, which allows to avoid scaleless integrals, must be carried out with 
the proper IR-subtractions. In the dimensional regularization such subtractions as a rule are performed 
automatically when one neglects the scaleless integrals, see for detailed discussion [53]. 
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2. Collinear and soft contributions in two-loop integral . 43) 



Now we return to the calculation of the overlapping soft and collinear contributions for the more 



complicate two-loop integral (43). From discussion in Sec III A it is clear that the situation with the 



soft-collinear overlap in (43 1 is quite similar to the considered above one-loop example. In the 2- loop case 
we also expect to obtain simple a logarithm InQ^/m^ originating from overlap of the collinear and soft 
regions. This can be written as 

+ UV-finite part 



c—p 



C—p 



Bt^ In 



(69) 



where Bj^^ denote appropriate part of the total coefficient _B''^ in (44). Computing B^^ we can follow 
the same line as in the one-loop case. 

We will perform all calculations in Z? = 4. In order to regularize the divergencies we introduce an 
infinitesimal gluon mass fj,. Such regularization looks quite natural in this case because the gluon mass 
plays the role of the virtuality cut-off for hard gluons in the hard subdiagram. Then the regularized soft 
contribution ( 46 ) can be written as 



) 03^*3 



dk 



W 6 



3J/33 



dk 



[kl 



A [kl 



(70) 



where we used that in Breit system p'_ ~ ~ Q, see Appendix A. Calculation of these integrals is a bit 
tedious but follows a basic line: two integrations performed by residue and the remnant integrals can be 
further computed keeping the most singular at — >■ terms. The details can be found in Appendix D. 
The result reads 



J. 



{2TTi) 



1 



TT 

"6" 



In 



(71) 
(72) 



where we again introduced r+ — ii^ lp'_ and r_ = 

Consider now the collinear term (|48|. We pick up from V only UV-finitc part (UV-f.p.) relevant for 
our calculation. Then the collinear integral reads 



UV-f.p. 



C(7_^ 



±^"3/33 



[-X3Q {k^ + k^ )] [-X3Q (fci + ^2 ) - ^2] [-xiQk^ - ^2] 



,(73) 



with 



*20(fc, ) 



[^[fhA^2fhMU I dk+dk+ dk,^dk2i_ 



— m?] [k\ — m?] 
2(p'fci) 



{p' -kl- k2f 



{p'-p',-k,r 



(k, + ka - p' + p',f {k, - p[r 



(74) 



The notation ^20 is introduced taking into account the structure of DA described in Eq.(|53[). Recall that 
according to ( 48 ) the integrals in Eq. ( 73 ) represent the convolution of the hard coefficient function with 
the perturbative DA ^'29. We carried out the calculations of the integrals in Eq.(74) keeping only the 
most singular contributions at the limit k~ ~ 0. The details are described in Appendix E. The result 
can be written as 



*2o(fcr) = [eii'uW2iVM] 



xe{0 < kY/Q < fcr™aJ^(0 < k^/Q < k^^^Jk^ In ( 1 + 

K2 



(75) 
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where dots denote the regular at the hmit ^ contributions which we do not consider for simplicity; 
the quantities ki2max ~ ^ denote the upper boundary for the relative light-cone components k~ /p'_. 



discussed in Sec. Ill A 



Their explicit value is irrelevant because we compute the integral in Eq. ( 73 1 with the leading logarithmic 
accuracy. Notice that from t he exp ression (75) one can see nonv anishi ng the boundary behavior of 'I'2o 

at /Cn 



Now we can substitute ( 75 1 into ( 73 ) and compute the regularized 
convolution integrals. For simplicity, we shall take into account the appropriate zero-bin subtractions 
changing the sign of soft contribution ([72|) in the sum ( 69 ) . Then we have 



-p ' 



where 



Jc—p' 



(27rz)2 


1 


1 




X 






1 


1 




X x^ 




(27ri)2 


1 


1 




X x^ 





dk. 



h ln(l + fci /k^ ) 

2 2 

'fcf-l-A:^) [k'^ + / xiQ\ 



dk^ 



dk'r 



. fc^'Mi + i/fc^') 

(1 + ^2^)' 



1 



TT 

"6 



In r_/Q 



The similar calculation for the second coUinear integral yields 

(27rz)2 1 



J c—p 



(76) 

(77) 

(78) 
(79) 

(80) 



Substituting obtained results into (|69[) and changing sign in front of the soft term we obtain 



s — c 



In 



{2mf 11' 



xixl yiyf 



TT 

IT 



In 



(81) 



We see that all regulators cancel as they should and we obtain a simple large logarithm. Recall that 
the coUinear contributions are associated with the hard rescattering term f[^\ Therefore this calcula- 
tion explicitly demonstrates that the soft and hard rescattering contributions are related and must be 
computed consistently. 

To summarize this section, we have demonstrated that the 2-loop diagram with massive quarks has 
a large logarithmic term originating from the overlap of the soft and coUinear regions. The appearance 
of this logarithm does not contradict to our factorization scheme. We have demonstrated that in the 
perturbation theory one can perform the consistent description of large-Q asymptotic of two-loop dia- 
grams using definitions of SCF (14) and DA ^ with free quarks instead of nucleons. We have also seen 
that coUinear integrals in the hard spectator contribution must be singular and such situation can be 
realized only with the help of additional assumption about the end-point behavior of the nuclcon DA. 



The factorization of the hard f[^^ and soft F'^"' contributions requires an additional regularization for the 
separation of coUinear and soft sectors. We carried out such separation for the case of 2-loop diagram but 
a realization of such scheme on a physical FF with nonperturbative DA and SCF is challenging because 
such matching involves different nonpertubative matrix elements. 



IV. DISCUSSION OF SCET FACTORIZATION SCHEME AND ITS 
PHENOMENOLOGICAL APPLICATIONS 



From the consideration of Sees. II and III we can conclude that the practical realization of the fac- 
torization scheme discussed in the introduction has some difUculties. Suppose that nucleon DA at low 
normalization point does not vanish at the end-points, as we discussed in Sec. III. Then the coUinear and 
soft convolution integrals are not defined in the both equations ([2| and ([s]). At a formal level we have 
established all the dominant regions and provided the definitions for the all nonperturbative quantities 
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but that is not enough. The comphcation arises due to the overlap of the coUinear and soft regions and, 
as a resuh, appearance of the specific "soft-coUinear" logarithms ^ InQ/A. Technically these logarithms 

(s) 

arise from the matching of divergent collinear and soft convolution integrals in the soft and hard 

F^'^'^ rescattering contributions. Therefore in order to describe corresponding logarithmic structure of the 
nucleon FF we must a formulate definite recipe for the separation of the collinear and soft modes. This is 
well known problem, for many processes where collinear factorization is broken due to singular end-point 
behavior of the convolution integrals see, e.g., [551 [Ml 

The presence of two hard scales: the hard ~ and hard-collinear ^ AQ in the description of nucleon 
FFs allows one to perform the factorization in two steps. First, one integrates over the hard fluctuations 
and passes from QCD to SCET-I. This effective theory includes hard-collinear modes which can be 
further factorized if the virtualities of the hard-collinear particles are quite large. Integrating over hard- 
collinear modes one passes from SCET-I to SCET-II which includes only collinear and soft particles with 
the virtualities of order A^. However if the value of Q is moderate, (for instance, the hard-collinear 
scale QA ~ ttiat ~ IGeV^ is not large in order to serve as expansion parameter) then one can not 
perform SCET-II factorization. Phenomenologically such situation is relevant for quite a large range of 
experimentally accessible values of momentum transfer. For instance, if is in the range of 4 — 25GeV'^ 
one obtains that QA varies between 0.8 — 2GeV^ assuming that A ~ 400MeV. Taking into account that 
QCD expansion of the jet functions is starting from we can expect that one needs much larger hard- 
collinear scale in order to see the dominance of this leading contribution. On the other hand we can not 

(s) 

ignore that the soft contribution F^ ' is suppressed by a Sudakov form factor. If this suppression is strong 
enough then then F^^^ would be subleading compared to the hard spectator contribution F^^\ However 
many phenomenological studies indicate that in the intermediate region of such suppression is still 
weak and the soft mechanism provides essential contribution to the physical FFs. Therefore in order to 
estimate the relative importance of the two terms it is necessary to include both of them consistently 
within SCET-I factorization scheme. But the end-point singularities in the hard scattering contribution 
make such program very complicated. One must introduce a factorization scale and factorize the end- 

(s) 

point contribution into F^^ . In some cases this difficulty can be solved using the universality of the 
SCET-I matrix elements. 

Suppose that we have different scattering processes which are described within SCET-I factorization and 
depend on the same SCET-I matrix element. Using universality of the SCET-I amplitude one can define 
the so-called physical subtraction scheme |47l I48| which allows one to perform the systematic calculations 
of the hard spectator scattering contributions associated with the symmetry breaking corrections. The 
idea of this approach is very simple: the SCET-I soft-overlap form factor or amplitude can be rewritten 
as a sum of one of the physical amplitudes and the corresponding hard spectator contribution. Then this 
combination can be used further for the analysis of physical amplitudes of other processes with the same 
SCET-I matrix elements. After such redefinition the end-point singularities in the combination of the 
hard spectator terms must cancel and one obtains the well defined hard correction. This scheme has been 
successfully used for analysis of different B-meson decay amplitudes and we expect that it can also be 
used for the analysis of the different hadronic reactions with the soft spectator scattering contributions. 

Let us illustrate the above discussion by one concrete example. Consider the following processes: 
7*iV — > iV describing proton and neutron form factors at large and 7*p — Tr^p describing wide-angle 
hard electroproduction of pion with s,f, ^ A^. We suppose that nucleon form factors are described 
by following tentative formulae: 

Ff (Q) = Ca{Q) { en /i"(Q) + /f (Q)} -f * * Hp * (82) 
F^Q) = Ca{Q) { en ftiQ) + ed /i"(Q)} + * * H„ * (83) 



where we used definitions (19 1,(20 1 and isotopic symmetry, symbols Hp „ denote the hard scattering 
kernel for the proton and nudeon cases, respectively. We also assume that the convolution integrals 
denoted by asterisk are regularized using some IR-regulator. Solving these equations with respect to 
SCET FFs one finds ( e„ = 2/3, = -1/3) 

A" = {2Ff + F^-<Sf* (2Hp + H„) * *} , (84) 
ft ^C^'{F? + 2^^r - * * (Hp + 2H„ ) * *} , (85) 



The pion production process can also be described as a sum of two contributions as shown in Figj5] So 
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+ 




FIG. 5: The sum of the soft and hard spectator contributions for the wide-angle hard pion electroproduction 



far we are not going to prove the exact factorization theorem for this process. Most important for us is 
that this configuration can be considered as one possible contribution to the nucleon helicity conserving 
amplitudes . The first term in Figjs] describes the soft overlap nucleon contribution and can be 
expressed in terms of the same SCET-I form factor The pion blob is described by pion distribution 
amplitudes. The second contribution in FigjS] can be associated with the hard spectator scattering and 
it is described as usually as a convolution of nucleon and pion DAs with the hard scattering kernel T^. 
Therefore schematic expression of the contribution in Fig[5] reads 



v47 = Lp^{z) * CH{z){eufi - edft) + *(a;^) * f^iz) * TH{x^,y^, z) * 



(86) 



where we again assume a some regularization for the divergent convolution integrals. Substitution of 
expressions for the SCET form factors (84) and (85) in Eq.(86) yields 



Af = Mz)*CH{z)/CA\{^Ff + AF^) 



(87) 



+ *(a;,) *(/3x(2) * [^Tff(a;„yi,z) - CH(z)/CA-{5Hp(xi,?/,) -t-4H„(x,,?/,)}j * ^(2/,). (88) 

Thus the soft overlap contribution is represented in terms of the physical FFs i^f'". The ratio Ch{z)/Ca 
depends only from the factorization scale associated with the evolution of pion DA. All Sudakov logarithms 
must cancel in this ratio. On the other hand, the end-point singularities in the hard scattering kernels 
must be compensated in the combination of Th and H„^p in rhs of ( 88 ) . The simple analysis show that 
the hard spectator correction is of order a^. At the same time the ratio Ch{z)/Ca = q;sC'lo(-z) +0{a1). 
The next-to-leading contribution can be computed from the one loop corrections to Ch{z) and Ca and 
the hard spectator corrections can appear only in the next-next-to-leading order. The more detailed 
analysis we are going to present in the separate publication. 

Therefore with this example we demonstrated how the SCET-I factorization on the soft and hard spec- 
tator contributions allows one to analyze the realistic hadronic processes at the intermediate momentum 
transfer. Of course, such method can be used when the number of observables is larger then the number 
of unknown SCET amplitudes. 

For the asymptotically large values of Q when the hard-coUinear scale is a good parameter for the 
asymptotic approximation one can try to perform the second factorization step and pass from the inter- 
mediate SCET-I to the low energy SCET-II. Performing this matching one must provide the solution of 
the soft-collinear mixing in order to treat correctly corresponding logarithms and avoid double counting 
between the soft and hard rescattering contributions. 

An interesting proposal for the separation of collinear and soft modes for the exclusive processes in the 
framework of SCET was discussed in |33| . In order to avoid double counting it is suggested to perform 
zero-been subtractions in the collinear convolution integrals. In our consideration in Sec. Ill we already 
used this idea in order to compute the investigated integrals in Z) = 4. The technique developed in [33] 
allows one to carry out zero-bin subtractions in dimensional regularization. In that case these subtractions 
remove collinear end-point singularities but generate the UV-poles. These UV-poles have been subtracted 
with the help of specific counterterms represented by non-local operators built from the collinear fields. 
Then, for instance, a divergent collinear convolution integral can be interpreted as following 



X0UO) 



0UO) lnp+/M 



+ 1 



(89) 
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In this example we used the pion DA (j)Tr{x) which, as it is often assumed, at the end point region behaves 
as (/)^(a; 0) ~ x. The auxihary scale arises after UV-renormalization and is quite similar to scales 
T± introduced by us in the calculations in Sec. III. 

Application of this scheme has been illustrated in [33j by a calculation of the transition form factor 
-Fp^. However, we expect that this analysis is not complete because the soft spectator contributions were 
not considered. The end-point singularities of the collinear integrals can be understood as a clear signal 
that such mechanism can also contribute in this case. We expect that, as in the examples considered 
in Sec. Ill, the soft spectator contributions will play the important role in the full compensation of the 
/i-i--scale dependence in the physical amplitudes. This means, that the logarithmic term in Eq.(89) 
with the derivative (/)^(0) is closely related to the soft rescattering contribution. Then one does not 
need the specific collinear UV-counterterms which were introduced in 33] . It means that the properly 
defined zero-bin subtractions must be unambiguously related with the UV-renormalization of the soft 
rescattering contribution. We expect that in this case one could carry out systematic analysis and avoid 
double counting. 

Formulation of such scheme is a difficult task and at present time we have not yet found a convincing 
technical realization of this idea. The most challenging problem is the consistent formulation of the zero- 
bin (or infrared ) subtractions for the collinear convolutions integrals with nonpertubative DAs (or other 
collinear matrix elements) and corresponding UV-renormalization of the soft rescattering contribution 
involving soft correlation functions, i.e. unambiguous matching of collinear and soft modes and conse- 
quently the definition of the hard and soft rescattering contributions. The other difficult moment is the 
absence of a well-defined regularization. Dimensional regularization can not be used for the calculation 
of the soft convolution integrals because corresponding integrands are scaleless. Therefore one has to 
invent the other method, for instance gluon mass, off-shell momenta or analytic regularization. As a 
rule, such regularization has potential problems beyond the leading order that makes difficult systematic 
calculations and analysis.'^ 



V. CONCLUSIONS 



In this paper we have investigated the soft spectator scattering term appearing in the calculation of the 
large-Q asymptotic form of the proton FF Fi. In our previous publication [1] we suggested that the full 
factorization formula for the nucleon FF Fi is given by the sum of two contributions: hard rescattering 
term F^^^^ and the second term involving soft spectator quarks which we refer as the soft rescattering (or 
soft spectator scattering) contribution f['^\ 

m • .... . . . fs) 

Using SCET technique in [T] we provided an explicit definition for the leading power contribution F-l 
in terms of the hard and hard-coUinear coefficient functions and the nonpertubative matrix elements. 
In this paper we have investigated more closely the proposed formula. We have computed the leading 
order jet functions J which appear in the matching from SCET-I to SCET-II. Using this result we 
have demonstrated that the corresponding soft convolution integrals in factorization formula ([s]) are 
logarithmically divergent in the ultraviolet and infrared regions simultaneously. This observation follows 
from the boost invariance of the soft correlation function S in Eq. ([3| and therefore does not depend on 
any model assumptions. It is natural to assume that such divergency arises due to the overlap of the 
collinear and soft regions. 

In order to clarify the situation we studied as example a two-loop QCD diagram with massive quarks 
and investigated the dominant contributions associated with the different regions. It was demonstrated 
that the overlap of soft and collinear sectors provides the large logarithm ~ InQ/m where m is the quark 
mass which plays the role of the soft scale. Therefore this confirms our assumption that the logarithmic 
divergence of the soft convolution integrals can be naturally explained by the overlap with the collinear 
region. 



During the preparation of this text new article |49| has been published where authors suggest a systematic formalism for 
the separation of soft and collinear regions and resummation of the arising large logarithms. In their terminology these 
are large rapidity logarithms. Application of this approach to the problems discussed in our presentation might be an 
interesting subject for future studies. 
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The obtained results allow us to make important conclusions about factorization FF Fi at large Q 
limit. We expect that the general scheme as described in Eqs ([2|-([3]) is valid but the separation of the 
soft F^'^^ and hard F^'^'^ rescattering contributions requires a certain prescription for the separation of 
the soft and coUinear degrees of freedom. From the overlap of the soft and collinear modes it follows 
that collinear convolution integrals in F^^"^ must be also divergent. This naturally leads to assumption 
that the nucleon distribution amplitude (fN Q-Q at some low normalization point has nonvanishing 
end-point behavior 

(90) 



(Pn{xi,X2 = 0, 1 - Xi) 7^ 0, (fN{xi = 0,X2, I-X2) 0. 



Such behavior is opposite to the asymptotic DA ipf^^Xi) ^ X1X2X3 vanishing at this limit. 

If Eq. ( 90 1 is realized then the corresponding collinear convolution integrals in factorization formula 



for Fj^'*' also suffer from the logarithmical end-point singularities which reflect the admixture of the soft 

mode. However, in the sum F^'^'^ and F^'^'^ the end-point singularities must cancel if computed consistently. 

The specifics of the factorization for Fi allows us to suggest that the correct treatment of the soft- 
collinear mixing becomes important only for a transition from SCET-I to SCET-II. In present paper we 
do not provide any systematic formalism for the factorization of the soft and collinear degrees of freedom. 
However, we would like to stress that the soft spectator contribution plays the important role in the 
correct description of such factorization. 

The picture with the soft spectator scattering mechanism naturally introduces a concept of two large 
scales: hard ~ and hard-collinear ^ AQ. This allows one to define the transition region or the region 
of moderate values of large momentum transfer Q. We suggest that it corresponds to the situation when 
the inverse power of hard scale l/Q^ is a good expansion parameter but hard-collinear scale is still not too 
large. Taking A ~ 300 - 400MeV and = 25GeV'^ one easily obtains that QA ~ 1.5 - 2GeV^. Hence 
in the situation when — 10 — 20GeV^ which overlaps with the majority of existed and upcoming 
experiments, one can perform consistently only the SCET-I factorization. Within such factorization 
scheme it is useful to take into account both soft and hard spectator scattering contributions because 
their relative contribution might be comparable if the suppression from the Sudakov form factor for the 
soft-overlap term is not sufhciently large. In order to perform a systematic consideration in this case we 
suggest to use the physical subtraction scheme which allows to express the SCET-I amplitudes in terms 
of physical ones. In such situation one can also to get rid of the end-point singularities from the hard 
spectator corrections. Applications to specific processes will be presented in subsequent publications. 
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Appendix A: Brief summary of used notations 



Through the paper we imply Breit frame 



/ 71 Tl\ 



and define the external momenta as 
,2 



^- Q=Q- 
Q 2' ^ ^2 




4m% 



(91) 



= Q + 0{m%/Q'), (92) 



where rriN is the nucleon mass. For the incoming and outgoing collinear quarks we always imply 



(93) 
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with the transverse momenta 



(95) 



and where Xi and x[ denote fractions of the corresponding momentum-component. Computing the 
Feynman diagrams in Sec. Ill we neglect by power suppressed components and assume 



P-Q^, P^-x^p, p ~Q-, p^~y,p, 



(96) 



In many formulas we use convenient notation Xi — I — Xi. We also use the following notation for scalar 
products 



(a • n) = a+ , (a • n) = a_ . 

and Dirac contractions 

Pfj.i'' = i> = p- 

Nucleon FFs are defined as the matrix elements of the e.m. current between the nucleon states: 



(pVeV(O) \p)=N{p') 

with nucleon spinors normalized as NN — 2mN 



7-(Fr + F.)-te^F. 



N{p), 



(97) 
(98) 

(99) 



Appendix B: Soft correlation function in perturbation theory 



The leading order perturbative expression for the reads: 



3 



) Lo = / '^^i *2 S{k+-uJi)d {k+ - Lo^) S (fcr - S {k^ - v^) 



iTr 

8 



k + 7_C (C7+ (k2 + m)) 



[tdlVl — k 



1± 



(100) 



(101) 



The factor 1/8 in front of trace is chosen for convenience. Calculation of the trace in the numerator 
yields: 



-Tr 



This allows us to write: 



(ki +rnjj-C (Cj+ (fca + 



-m^ + (fcij_ • k2±) 



3m2 1 



LO 



167r8 4 



dki± dk2±_ 



[ujiiyi — fcjj_ — + is] [uj2V2 — fc|_L ^ + 



(102) 



(103) 



In order to proceed further we must take into account the specific properties of the jet functions. We 
always assume that the soft fractions uJi and Ui are positive. Mathematically it comes out from the 
analytical properties of the diagrams and imposes specific restrictions on the integrand in Eq.(103|. In 



order to see this assume that — oo < Wi < oo but then we keep the Feynman ie-prescription in the jet 
functions. From calculations of the diagrams in Fig[2]one can easily obtain that all denominators of the 
jet functions in Eqs. (29)-(32| are defined with —ie, for instance 



1 



1 1 



1 



1 



1 



ll>i + 0J2 L^i ll>2 1^1 + ^2 ~ *e] ['^i ~ is] 1^2 — is] 
The same arguments also true for i^i. Consider now the convolution integrals 

1 1 1 



djjjiduj2-. 



cji + UJ2 — *£ Y^l — — i£\ 



(104) 



(105) 
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where Svii^i) is represented by expression (103 1. Computing diOi by residues we obtain that nontrivial 
results originates only from the poles of the propagators in S'y(cji) in (103). Therefore this allows to us 
to represent the propagators in ( 103 ) as (5-functions: 



[ujii>i — fc^j_ — m? + ie][uj2'^2 ^ k^j^ — m? + ie 



-^{2mfe{uj, > 0)9{iy, > 0) 



Then we obtain 

( S^y') 



_ 3m2 
LO IBtt^' 
3m2 



167r6 



(106) 



X S (wii^i - klj_ - w?) 6 {uj2V2 - - w?) . (107) 



{v, > 0)e{uj, > 0) e{uji?]i > m^)0{uj2r]2 > m^)- (109) 



Appendix C: Derivation of the collinear contribution D 



c — p' 



Consider first denominator ( 43 1 . In the collinear region ( 47 1 we obtain: 



[{p-ki-ka) - m^] {p ~ p3 - ki ~ k2) (p-ps-^i) [(pi - ^i) - m^] 

^ [~Q {K + ^2)] [-Qx3 {K + ^2)] [-QxzK] [-Qxik^] 

c^Q^{xlx,)[- {k^+k^)]'[-k^]\ 
The remaining propagators are soft, of order A'^. Hence for the denominator we obtain 

Den ^ Q«Ai2 

In the numerator we have 

Num = eU^i + [7'' ei] ^"2 1\P' ~P3-ki+ m)7"(fc2 + m) {p - Pi - hW 6 



ki ^ k2 m 



7np-fci-fc2 + m)7^ e 



where we single out by brackets [...] the numerators of the hard propagators and hard gluon 
Using that ~ we can rewrite this piece as: 

[7^6]® 1^ ip - Pz - kiW ^2 ® j1 (p ~ ki - k2 + rnj 7'^ ^3 

c.(-4)(pA;i)[7^ei]®[7^'6]®W6]- 

Then we obtain 

Num = 2Q{~kY)^W {ki + m) [7^"^] e^7'(p' - p'3 - k + m)rCk2 + m) [j^ ^2] 
xfW"(p'-fci-fc2+rn) [7'^ 6]. 

Recall that ~ Q, then for the remaining terms one finds 

CiY {ki±) [7^a] C3 r - hi) bi 6] C2 i'i-ki±)rik2±) 6] 

^A4^"iriei?^r26elr36, 

where Ti denote certain Dirac matrices. Hence one obtains: 

Num - Q^A^, 



(110) 
(111) 

(112) 

(113) 

(114) 
(115) 

vertices. 



(116) 



(117) 



(118) 
(119) 
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and for the whole diagram one finds: 



The coUinear contribution can be written as 



X 



[-Q{k^+k^)]'[-Qk^] 



with 



(120) 



(121) 



dk^dk^ dki±dk2±-r-^ 57770 ^ 

{^'27* [(yap' - fci) + ™] 7"(fc2 - m)}^ 7" [(p' - fci - fca) + m] }^ 



ip' - fci - ^2)^ - m2 (fci + fc2 - yaP')^ (j/aP' - h f - rn? 



(fci - yip') 



(122) 



Eq.(122| yields the c ontri bution to the evolution kernel at 2-loop approximation. The integral with 
respect to k^2 i^i Eq-(121| can be interpreted as a convolution integral of leading order hard coefficient 
function with the given part of evolution kernel. 



Appendix D: Calculation of Ds 



The expression for the soft integral in Eq. ( 72 1 reads 



Js — -2 -2-^-5' (123) 



with 



= y dk2±dki± J dk2 dk^ ^ 



k+ - T+/y,] [-{k+ + k+) - T+/y^] [-{k+ + fc+)] 



(124) 



1 Tfl^ 

'^^^ '^^^ [-fcr - r_/a;i] [-(fc^ + k^) - r^/x^] [-(fcf + ^2")] [fci - m^] [fc? - m^] ^^^^^ 
Let us redefine the notation as 

fc+ Pi, kl _> ai, dki± -> dki, ti = t^/xi, T3 = t^/x^, t[ = r+/yi, Tg = r+fy^. (126) 
So that 

1 



Is= J dk2dki I dPi dP2 
X / daida2 



[-Pi - n] [-(/?! + f32) - r-s] [-(/3i + ^2)] 
1 



[-ai - r{] [-(ai + 0:2) - T3] [-(ai + £12)] [a2;92 - - m^] [ai/3i - kj - m?] ' 

Expressions in square brackets implies [...] = [... + le]. Next we use the same trick as in Appendix A: we 
integrate over dai 2 by residues, rewrite the poles as (5-functions: 

^ {2Tiife{Pi)e{ai)5{ail3i - kj ~ 771^)6(02^2 - kl - m^) (127) 



[q!2/32 — ^2 — fn^] [aiPi — k\ — m? 

and integrate over transverse momenta. This yields 

1*00 

Is = {2nifm^ / dPi dP 
Jo 

J 

daida2 



[Pi+n] [(13^ + 132) + T3] (/3i+/32) 
0{a2l32 > m^)9{aiPi > m^) 



[ai + t[] [(«! + 02) + T3] («! + 02) 



23 



After simple substitutions this integral can be written as 



Is = {2mfm^ / 
Jo 



dl32 



1 



daida2 



' " (1 + /32) 7o " ("1 +^2) 

^ r , . Pi 1 1 



Using that 



/3i 



T3 



[/3i(l + /32)+r3] (I + /32) (l + /32)[/3i(l + /32) + r3] (l + /32)^ 



(128) 



1 



1 



(ai + /3ir() [ai + q;2 + /SiTg] (ai+a2) 

1 



1 



(ai+/3ir{) [ai + a2 + /3iT;^] 
1 



(ai +^2) (ai +/3it{)' 

where we neglected small contributions proportional to infinitesimal mass. Therefore we obtain 

1 



(129) 
(130) 



Is = {2TTiYm^ / d/3; 



g(ai>m2) g(a2/?2>m^) 1 ^^^^^ 

-(l + /?2)^7„ '"^'"^ (ai+a2)^ io (/3i+r(n) (o,+A)^'''^ 



(2.^)^m2 /°° da. / da 



1 



«! — r^Ti 7o (ai + 0:2)^(1 + ?™^/a2) 

(2..)^ r ^ ( mar +ln K^n]) ^^^-^^ 
Ji ai (1-ai)^ 



Therefore we finally arrive at: 



Is = {2nif In [r+r_/m2] ( 1 _ _ ) + o{l). 



(132) 
(133) 
(134) 

(135) 



Appendix E: Calculation of ^2o{k^ 



Let us rewrite the Eq. ( 74 ) as 



*20(fc-)=Kiy]a,Ki7']aJe3]o3-/20, 



with 



J20 = j dk^dk^ dki±dk2± 



m 



-2(p'fci 



(p' - fci - fc2)^ - (p' - P3 - fcl)^ - TO^ (fc^ + fc2 - P' + pjj)^ (fci - P'l)^ 



(136) 



(137) 



The exact answer for J20 is very complicated. We shall compute this expression only in the region of 
small k~ — >■ assuming that their ratio is fixed k^ /k2 ^ C(l)- Again, redefine for simplicity light-cone 
decomposition: 



= K 2 + 2 + ^«-L = 2 ^'2 



(138) 
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and rewrite Eq.(|137|) as 



J20 = / dPi df32 dk2dki 



QPi) 



1 



(139) 



[(A + P2) (ai +a2-Q)- (kl + /C2)2] + 132) («! + ^2 - ysQ) - (^i + h?] 

Performing integrations over /3i and (32 by residues we can neglect by all poles for which 

at > Q - aj or > Qyt - aj, (140) 

because in this case can not be small. Then we have contribution only from the poles associated with 
the first two propagators in Eq.(139) (in the first line). The result can be written as 

J20 ^ {2mf0{O < ai < QaT") ^(0 < ^2 < QaTl 

X j dPidl32dk2dki5{a2P2 - k^ - m'^)S{ail3i ~ kj - m^) 
with the denominator 



(-Q/?i) 
Den 



(141) 



Den = [(ai - yaQ)/?! - kl] [(ai - yiQ)pi - kl] [(A + P2) (ai + 02 - Q) - {ki + k2f] 

X [(/3i + /32) (ai + «2 - ^sQ) - (fci + k2Y] (142) 

The maximal values a^^^ ^ a™^'' ^ 1. Their explicit values are not important. The integrand can be 
further simplified. We can neglectthe small a.i <^ Q in the denominator (142): 

Den ~ [(-ysQ)/?! - kl] [i~yiQ)pi ~ kl] 

X [(/3i + M i-Q) - (fci + k2f] [(/3i + P2) i-mQ) - {kl + k2f] . (143) 

Then we take into account that the dominant contribution arises from the region where 

kf = a^ft - < yiQ, (144) 



as it follows from the ^-functions in Eq.(141). Therefore we can also neglect the transverse momenta fc^ 
in the propagators (143 1: 



(— m^) 



1 



Den - [(-y3Q)A] [{~yiQ) l^i] , , ^ , riM r^fl 1 - " -2^3 0^0 o ^2■ 



(145) 



Finally we obtain 

J20 ^ {2Tiif0{Q <ai< gar") ^(0 < "2 < Qar") 
X y d/3i dP2 dk2dki S{a2(3; 
A simple calculation yields: 



(-1) 



I - TO^) (5(ai/3i - A? - m^) 



^1 (/3i + /32) 



(146) 



J20 ~ (27^^)26'(0 < ai < QaT"") 6'(0 < 0^2 < Q^r'') 



(-1) 

2(03 



a2 ln(l + ai/a2) 



^ {2mfei0 < k^/Q < fc^„_) 6(0 < k^/Q < fc^^J 



k^Hl + k^/k^). (147) 



Substituting this into Eq.(136) we obtain the required result. 
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